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Abstract. Polymer quantum systems are mechanical models quantized similarly as loop 
quantum gravity. It is actually in quantizing gravity that the polymer term holds proper 
as the quantum geometry excitations yield a reminiscent of a polymer material. In such an 
approach both non-singular cosmological models and a microscopic basis for the entropy of 
some black holes have arisen. Also important physical questions for these systems involve 
thermodynamics. With this motivation, in this work, we study the statistical thermody- 
namics of two one dimensional polymer quantum systems: an ensemble of oscillators that 
describe a solid and a bunch of non-interacting particles in a box, which thus form an ideal 
gas. We first study the spectra of these polymer systems. It turns out useful for the analysis 
to consider the length scale required by the quantization and which we shall refer to as poly- 
mer length. The dynamics of the polymer oscillator can be given the form of that for the 
standard quantum pendulum. Depending on the dominance of the polymer length we can 
distinguish two regimes: vibrational and rotational. The first occur for small polymer length 
and here the standard oscillator in Schrodinger quantization is recovered at leading order. 
The second one, for large polymer length, features dominant polymer effects. In the case of 
the polymer particles in the box, a bounded and oscillating spectrum that presents a band 
structure and a Brillouin zone is found. The thermodynamical quantities calculated with 
these spectra have corrections with respect to standard ones and they depend on the poly- 
mer length. When the polymer length is small such corrections resemble those coming from 
the phenomenological generalized uncertainty relation approach based on the idea of the 
existence of a minimal length. For generic polymer length, thermodynamics of both systems 
present an anomalous peak in their heat capacity Gy. In the case of the polymer oscillators 
this peak separates the vibrational and rotational regimes, while in the ideal polymer gas it 
reflects the band structure which allows the existence of negative temperatures. 
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1 Introduction 

In coping with the challenge of quantizing gravity the loop quantization approach [1, 2] has 
proved convenient to incorporate the background independent character demanded by general 
relativity. Important progress in this approach include the avoidance of the classical singularity 
which in loop quantum cosmology is replaced by a quantum bounce [3] in physically motivated 
models, [4, 5, 6] , and a microscopic basis for the entropy of some black holes that is in accordance 
with the Bekensteim-Hawking's semiclassical formula [7, 8, 9]. 

However some further physical questions in regard to cosmology and black holes necessarily in- 
volve thermodynamics. As it is well known primordial particle backgrounds including neutrinos, 
gravitons and photons, originating in the early universe provide windows to explore such early 
stage [10, 11]. For instance the stochastic graviton remnant and its statistical properties have 
been studied [12, 13]. Also the different contributions to the spectrum of gravitons produced 
during the super-inflationary period have been investigated [14, 15, 16, 17]. In regard to black 
holes, from the loop quantum gravity perspective, the challenge remains of describing black hole 
evaporation including in particular its thermodynamical aspects (see e.g. [18] for recent work.) 

Rather than dealing with the thermodynamics of loop quantized gravitational systems a more 
tractable problem is to consider the statistical thermodynamics of polymer quantum systems 
[19, 20, 21]. The latter are mechanical systems quantized following loop quantum gravity. Here 
an important comment in regard to the polymer term is in order. In the gravitational case, 
loops, or more strictly, graphs, label quantum states of gravity. This yields a picture resembling 
polymer materials, which justifies adopting the term polymer quantization for gravity. However, 
as we will see below, for mechanical systems states will be labeled by point sets belonging to 
a lattice. Thus, although the term polymer is inherited from the gravity case it is not actually 
realized in the mechanical one. We should stress at this point also that in loop or polymer 
quantization a length scale is required for its construction, while for the gravitational case this 
is identified with Planck's length, in the mechanical case it is just a free parameter and we refer to 
it as the polymer length scale. It should be mentioned though that, technically, the quantization 
here dubbed polymer, was previously considered from a different perspective in the form of a non 
regular representation of the canonical commutation relation [22]. Interestingly, a quantization 
based on difference operators has also been considered [23]. 

Polymer quantum systems have been convenient to illustrate some features arising in loop 
quantum gravity [19, 24]. In particular they have the same configuration space as that of loop 
quantum cosmology [25]. The continuum limit of the polymer quantum system has also been 
explored using ad hoc renormalization schemes in which the polymer length scale runs [20, 21]. 
Investigating whether polymer quantum system admits Galilean symmetry has been reported 
in [26]. Inspired by the cosmic singularity avoidance this quantization has been used to explore 
potentials such as 1/r [27], and 1/r 2 [28]. 

In this work we shall study the thermostatistics of two simple polymer systems, namely, an 
ensemble of oscillators and a bunch of noninteracting particles in a box. The paper is structured 
as follows. In Section 2 we review the basics of polymer quantization of mechanical systems, 
in particular the eigenvalue problem. For the harmonic oscillator the corresponding eigenvalue 
problem can be casted in Fourier space as a second order differential equation. In this way 
one can see its spectrum is just that of the standard quantum pendulum. Two regimes can 
be seen to appear: an oscillatory or vibrational one and a rotational one. In the first, the 
oscillator in Schrodinger quantization is recovered, at leading order, while in the second regime, 
the polymer effects are dominant. Furthermore, we solve exactly the eigenvalue problem for the 
polymer particle in a box that takes the form of a second order difference equation. We obtain 
a bounded and oscillating spectrum that features a band structure and a Brillouin zone. 

In Section 3 we calculate the corresponding thermodynamical quantities with these spectra. 
They depend on the polymer length and behave differently from those in standard thermo- 
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dynamics, which is recovered when the polymer length is considered small. In this case the 
thermodynamical variables can be written as a series on the small polymer length, similarly 
to what happens in the phenomenological Generalized Uncertainty Principle (GUP) approach 
based on the idea of the existence of a minimal length [29, 30]. These were applied previously to 
an ideal gas [31, 32, 33], and radiation [34, 35]. (The modification to the statistical mechanics 
of systems were also studied from the perspective of the extension to the Standard Model that 
have Lorentz violating terms [36], and the case of radiation was also studied with corrections 
arising from loop quantum gravity [37, 38].) In some sense this behavior could be expected since 
the uncertainty relation for polymer systems is quite similar to GUP [19, 39]. We shall show 
quadratic corrections occur in the energy of the polymer particle in a box, just as in GUP, but 
with opposite sign. For a generic, not necessarily small, polymer length, thermodynamics of 
both systems present an anomalous peak in their heat capacity Cy . For the polymer oscillators 
this peak separates the vibrational and rotational regimes, while in the ideal polymer gas it 
reflects the band structure. It is worth stressing that, for the gas, the band structure allows for 
the existence of negative temperatures. 

Finally in Section 4 we discuss our results and point out some perspectives. 

2 Eigenvalue problem in the polymer representation 
of quantum mechanics 

In this section we describe the main features of the polymer quantization of a non relativistic 
particle moving on the real line [19, 20], and describe briefly the corresponding eigenvalue 
problem. To do so we start by noticing that instead of the Heisenberg algebra, involving posi- 
tion q and momentum p of the particle 

[q, q] = [p, P] = 0, [q, p] = qp - pq = il, 

where I is the identity on Hilbert space T~L, one adopts the Weyl algebra 

tf(Ai) • U(X 2 ) = U(\i + A 2 ), V(vn) ■ UM = V(jjli + fj, 2 ), 

According to Stone-von Neumann theorem the elements of the above algebras can be related 
under certain conditions, including in particular weak continuity. Namely 

which however, does not hold in the polymer case. 

In the usual or Schrddinger representation of quantum mechanics the Hilbert space is T-L = 
L 2 (M.,dx) with the Lebesgue measure dx. Instead, in the loop or polymer representation the 
kinematical Hilbert space % po iy is the Cauchy completion of the set of linear combination of 
some basis states {(xj)}, whose coefficients have a suitable fall-off [19], and with the following 
inner product 

1 f a+T 

(xAxj) = lim — / dke ik ( Xi ~ x ^ = 6 X . x ., 

where S XijX . is the Kronecker delta, instead of Dirac delta as in Schrodinger representation, then 
we say that the orthonormal basis is discrete. The kinematical Hilbert space can be written as 
Hp \y = L 2 (Mrf,d/id) with d/id the corresponding Haar measure, and the real line endowed 
with the discrete topology 1 . 

1 In the momentum representation, the configuration space is the Bohr compactification of the real line Rs, 
[19, 25]. 
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As we already notice Stone-von Neumann's theorem is no longer applicable since the opera- 
tor V(fi) fails to be weakly continuous on the parameter \i due to the discrete structure as- 
signed to space [40, 19]. The shift operator V[p) is not related to any Hermitian operator as 
infinitesimal generator. Hence, for the representation of the Weyl algebra we choose the posi- 
tion operator x and the translation V(/i) instead of the momentum operator. Its action on the 
basis is: 

x\Xj) = Xj\Xj), V(fJ,)\Xj) = \Xj — /i), 

fulfilling also: 

[x,V(n)] = -»V(n). 

Since there is no well defined momentum operator, any function on phase space which depends 
on the momentum has to be regularized. In particular this is so for the Hamiltonian. To do so 
we introduce an extra structure, namely a regular lattice with spacing length fiQ. This analogue 
of what happens in Loop Quantum Cosmology, where there is a fundamental minimum area 
[3, 4, 5, 6] given in terms of Planck length. 

Let us consider /io > as any fixed scale. In general /io can be function of x but here 
we suppose it is constant. One of the simplest options to define an operator analogous to the 
momentum operator is: 

K m = ^-{V{^)-V(-^)). (2.1) 

This choice can be thought of as the formal replacement 

P -> — sin (/Hop), 
Mo 

where the right hand side is given by (2.1). Some authors have studied a semiclassical regime in 
which the expectation value of the Hamiltonian is taken with respect to a semiclassical state to 
yield an effective dynamics that can be seen, at leading order, as obtained from the replacement 
p — > sin (nop) I no in the classical Hamiltonian. This has been proved useful in several models 
[41, 42, 43, 44, 45]. 

Thus the polymer Hamiltonian is written as 

H w = ^2 t 2 " - n-Mo)] + W{x), (2.2) 

where H^(a;) is a potential term. The dynamics generated by (2.2) decomposes the poly- 
mer Hilbert space T-L po \y, into an infinite superselected finite-dimensional subspaces, each with 
support on a regular lattice 7 = 7(^0, xo) with the same space between points where 
7(/xo, xq) = {nfJ-o+xo I n e Z}, and xq £ [0, /io). Thus, choosing xq fixes the superselected sector. 

There are at least two possible ways to regain Schrodinger quantum mechanics. One is 
to consider the polymer length scale ^0 as small such that the difference equation (2.2) can 
be approximated by the usual differential Schrodinger equation [19]. Another possibility is to 
consider the introduction of the lattice as an intermediate step, after which it is necessary to 
carry out a renormalization procedure, as it is usually done in lattice theories [46]. In polymer 
quantization the final result turns out to be the usual Schrodinger quantum mechanics [21, 20]. 

Next we consider the eigenvalue problem: 

H ltQ \ij;)=E\^). (2.3) 
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In the case of the Hamiltonian (2.2) on the lattice j(xo, Mo)> an Y state \ip) G % po i y is of the form: 
W = ^2^X0 + 3 Ho)\x + 3 Ho) ■ (2.4) 

The substitution of (2.4) into (2.3) gives a difference equation in the position representation 

h 2 



2 [Win 3) ~ ^( x j + Mo) - - Mo)] = - VF(zj)] ^(xj). 



2m^ 

Using the Fourier transform defined by 

fP(k) = (£#) = J2^o+3^)e- ik ^ 0+j ^\ 

where k G [— tt/^o, tt/mo] and with ^(fc) satisfying the condition 
VMo/ 



e M o tp 

' Mo 



the equation (2.5) reads: 



h 2 



E — cos kfXQ ip(k) 



in 



h 2 



-ifc(x +j> ) 



3& 



(2.5) 



(2.6) 



In our analysis two specific cases will be considered: the harmonic oscillator and the particle in 
a box. 



2.1 Harmonic oscillator 

The polymer harmonic oscillator has been already studied in [19] and [47] 2 , using the potential 
W(x) = huix 2 /2d 2 , with d 2 = h/mu a characteristic length of the oscillator. After performing 
the Fourier transform in the r.h.s. of (2.6) and using the quadratic potential, one obtains a second 
order differential equation which can be recognized as a Mathieu equation [49] 



+ (a - 2< 7 cos20)V>(</>) = 0, 



(2.7) 



where 



kflQ+TT 



= j5 — lj and q = 4A , and we introduce A := fio/d as a dimen- 

sionless length parameter. This system is just a quantum pendulum in k space [20]. For low 
energies one recovers the harmonic oscillator behavior, while for high energies it becomes a free 
rigid rotor [50]. Equation (2.7) has periodic solutions for particular values of = a n ,b n called 
the Mathieu characteristic functions, that depend on n, [49]. The corresponding wave functions 
can be written in terms of the Mathieu elliptic sine and cosine [47] . The energy eigenvalues can 
be expressed as follows [49, 47] 



E% n — 



flUJ 
J 2 



1 A 4 [A 

1 H On 77 

8 \A 4 



>7 



A 



1 + — b n+1 



( 4 



\A 4 



(2.8) 
(2.9) 



2 Actually a previous study appears in [48], in a different context, transport of electrons in semiconductors. 
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An asymptotic expansion for the characteristic functions, assuming //q <C d yields that the 
energy spectrum E2 n ~ ^n+i, can be approximated as [19, 47] 

E n = (n + I) tiu - ( 2n2 \l n + 1 ) X'hoo + 0(A 4 ). (2.10) 



2 J V 32 

As expected the spectrum consists of the standard harmonic oscillator plus corrections of or- 
der 0(A 2 ). However, this spectrum is not bounded from below. To enforce correspondence 
with the Schrodinger quantization, we can ask that the leading term be larger than the first 
correction. This gives us a maximum n max that depends on A, 

n max « A" 2 . (2.11) 

This means that, for small A, we can not probe the spectrum with values of n greater than those 
allowed by (2.11). As in [19] with the values of /j,q = lCP 19 m, corresponding to the maximum 
experimental attainable energy today, and with d = 10 _12 m for the carbon monoxide molecule, 
one gets A = 10 -7 . In this case n max ~ 10 14 . This is consistent with [21] and [20] where a cut-off 
in the energy eigenvalues was introduced that depends on the regulator scale in order to perform 
the renormalization procedure that is necessary to implement a continuum limit of the theory. 
On the opposite limit for A S> 1 the eigenstates are Em ~ -E^n-i, for n = 1, 2, . . . 

£ n = ^ + ^n 2 + 0(A- 6 ), (2.12) 

where the ground state depends only on A -2 , and actually falls off as A increases. As pointed 
out in [47] this case may be relevant for the cosmological constant problem. For the excited 
states ii / 0, the A 2 term dominates. In this regimen it is better to interpret A as a ratio of 
energies. Indeed the dimensionless parameter A 2 = , where E osc = fooo and E po \ y = ^2 

that corresponds to the coefficient of the kinetic term in (2.5). Then this regime is the case 
of E osc S> E po \ y ; if we consider the polymer length to be the Planck length then this case 
corresponds to the trans-planckian region relevant for both inflationary cosmology and black 
holes [51]. 

Following the comparison with the quantum rigid rotor [52], we recall that the spectrum of 
a quantum rigid rotor in two dimensions is ^T - , with / = being the moment of inertia 

and R the radius of the rotor. Notice that the second term in (2.12) has the same functional 
dependence on n as for the rotor with effective moment of inertia J e ff = an d -R e ff = >~d. 

2.2 Particle in a box 



The next example is a particle confined in a box of size L = NfiQ. The free polymer particle 
was first studied in [20]. Here however we confine it to a box. In this case, instead of working 
in Fourier space we can use directly the difference equation (2.5). As in the standard case the 
potential is defined as 

W{ Xj ) = {^ ^< Xj < XQ + L, (2i3) 
I 00, otherwise, 

where xq £ [0, fj,o). This potential is then realized through appropriate boundary conditions over 
elements of H po iy The particle behaves as a free particle inside the box and vanishes outside, 
namely 



?/>(x ) = ip(L + xo) = 0, Vi 6 [0, Ho)- 



(2.14) 
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We can conveniently rewrite equation (2.5) by using Xj = xq + jfj-o, as 

2mE$ 



<KJ + 2) 



h 2 



V>(j + 1) + i/>(J) = o. 



(2.15) 



Following [53], we propose the solution of the difference equation of second order (2.15), to be 
= a\r{ + a 2 r 1 2 , (2.16) 
where a% are constants coefficients and rj are the roots of the characteristic equation: 



r2 _ {2 _2rnE^ , r + 1 = a 



h 2 



The solutions of (2.17) are 



(2.17) 



1 ± 1 / 8m£/x 2 / mEy 2 ^ 



h 2 



2 V fr 2 V 2/i 2 



For positive energies, the argument of the square root gives us a relation between energy and 
2h 2 /m^iQ. If E > 2/i 2 /m//Q, then the roots r± are real numbers 3 , but incompatible with the 
boundary conditions, therefore they yield the trivial solution ip(xj) = 0. This leads to the idea 
that the minimum length scale /io imposes a cut-off on the energy, and hence, energies greater 
than the cut-off are unphysical. Thus, the only meaningful physical case is when E < 2h 2 /mfj,Q 
which gives us complex roots for r±. The solution is given by: 



il>(j) = Ci S m(j6) + C 2 co S (j9), 
which is a parametrization of (2.16) in polar coordinates, 



(2.18) 



COS I 



1 



mEfiQ 



sin 



1 8mEfj, 2 ( mEji 2 \ 
Q -2\—^\-^- X )- 



Imposing the boundary conditions (2.14) in (2.18) we find that 6 = nnX. In this case A 
Ho/L = 1/N and n G Z. The eigenfunction of (2.15) turns out to be 



ifinfaj) = C sin I mrX 



Mo 



C sin ( nir 



3_ 
N 



0<j <N, 



(2.19) 



where C is a normalization factor 

N 



c 



^sin 2 (rnrXj) 

3=0 



With (2.19) we can write the eigenstate (2.4) as 
\ipn) ■ \/ — > mii I n~A— ! \.r,} 



\ y^sin ( nirX— ) 



(2.20) 



3 One of these cases give rise to degenerate roots, such that the solution (2.16) changes for ip(j) — r J (ai +a.2j) 



8 



G. Chacon-Acosta, E. Manrique, L. Dagdug and H.A. Morales- Tecotl 




Figure 1. The solid line (blue) corresponds to the first Brillouin zone for the spectrum (2.21) and 
the dashed line (red) is the second order approximation equation (2.22). The dot-dashed line (green) is 
the case reported in the GUP literature corresponding to equation (2.23) which presents and opposite 
tendency with respect to the previous one. The dotted line (black) corresponds to the standard case 
A = 0. 

We notice that the sum in (2.20) goes form j = 0, . . . , N. Since n is an integer it is clear that 
we can not build the n + 1 eigenstate because it would depend on the previous n states. Thus 
< n < N. The corresponding energy spectrum is found to be bounded [54] 



ri 2 



(1-cosnvrA), n 6 {1, 2, . . . , N - 1}. 



(2.21) 



The energy spectrum (2.21) resembles the tight binding model for particles in a periodic potential 
that is not infinitely high and allow tunneling with the nearest neighbors sites [55, 56, 57]: 

E(k) = E* — 2 A cos k/j,q, 



where E* = and A represents the relevant non diagonal terms that give rise to the energy 

band, n is interpreted as a wave vector which take values on [— ^], and which, by choosing 
periodic boundary conditions [56, 57], it takes discrete values k = £p with N G Z and —N < 
m < N. 

Hence, the polymer particle in a box is analogous to the tight binding model of a particle 
in a periodic potential with periodic boundary conditions, having band energy A = E*/2. So 
there is an energy band appearing in this polymer system given by a energy spectrum bounded 
from above and below, and a kind of Brillouin zone as shown in Fig. 1. 

If we expand cosnvrA, for A< 1, i.e. large size of the box as compared to the lattice spacing, 
up to second order, we get 



E n 



h 2 n 2 TT 2 



h 2 n 4 ir\,. 



2mL 2 24mL 2 



(2.22) 



Using fiQ = 10~ 19 m and with the approximated spectrum (2.22), the correction will be significant 
only when n ~ 10 17 . 

The energy spectrum of a quantum particle in a box has also been studied in the framework 
of a GUP [58]. In this frame the modifications induced by a minimum length £ m [ n on the wave 
function and the energy spectrum of a particle in a one- dimensional box, result in a modification 
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of the Schrodinger equation which transforms it in a fourth order differential equation. The 
energy spectrum has a correction proportional to the squared of the minimum length 01 



mm 



n 2 ir 2 h 2 i 2 - n 4 7r 4 h 2 

E n ' = IT— FT + -ff o T2 - 2 - 23 
2mL 2 L 2 3mL 2 

This result can be compared with (2.22). We can see that the dependence on the minimum scale 
is quadratic in both cases but they feature a sign difference thus signaling opposite tendencies. 
As we already mentioned this coincidence is not surprising since polymer systems have similar 
modifications to GUP in their corresponding uncertainty relation [39]. However, in [58] the usual 
boundary condition for a particle in a box is used; here, since the measurement of the position 
has some uncertainty proportional to the minimum length, the position of the walls is not 
determined precisely. This problem is solved recalling that in polymer quantum mechanics the 
dynamics is defined on superselection spaces for which there is, in each one, a proper boundary 
conditions defined by (2.13) and (2.14). 

3 Polymer corrections to thermodynamic quantities 
of simple systems 

In standard statistical mechanics [59], the canonical partition function Z is defined as the sum 
over all possible states 

Z(/3) = ^exp(-/3£ n ), (3.1) 

n 

where j3 = (kT)^ 1 , k is the Boltzmann's constant, and T is the temperature. With it we can 
calculate all thermodynamical quantities of the system, through the definition of the Helmholtz 
free energy [59] 

F = -jlnZ, (3.2) 

where TV is the particle number. The relationship with other thermodynamic quantities such as 
the equation of state, entropy and chemical potential can be obtained by the standard relations 

dF , . 

"=-8L- < 33 > 
dF 

Moreover the internal energy and the heat capacity can also be related to the partition function 
as relations 

Kr d\nZ MdZ . . 

u= -"-ar = -zw (M) 

BTT 

Cv = -W 2 -Qp- (3-7) 

Thus, all we need is to calculate the partition function (3.1) using the corresponding energy 
spectrum. 

Here we use both closed forms (2.8), (2.9) and (2.21), as well as the approximate forms 
(2.10) and (2.22) of the spectra to calculate the canonical partition function (3.1) and to obtain 
thermodynamical quantities for the polymer solid and ideal gas. 
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In this work we have adopted the Maxwell-Boltzmann statistics for particles with no spin. 
However, it has been argued that modified statistics may be needed in quantum gravity and 
discrete theories. The relation between the discrete structure and statistics is still open [60, 61]. 



3.1 Ensemble of polymer oscillators 
3.1.1 Exact spectrum 

To calculate the partition function using the spectrum (2.8), (2.9), we split the sum into two 
parts 



Z((3) = ^exp 



(3kuj ( A 4 
' 1 H a 



A 2 



8 " V A 4 I 



cxp 



Pfuj ( A 4 



(3.8) 



the first term takes into account the even, and the second term the odd, parts of the spectrum 
with n and n' running over all integers. 

To determine the thermodynamical quantities we first write the Helmoltz free energy as 



F 



ln i e ^E 



A 2 A / A 2 

exp I —[3huj—a n I + exp I —(3hoj—b n +i 



where we omit the argument of the characteristic Mathieu functions to avoid cumbersome 
expressions. From (3.3) and (3.4) it can be noticed that the equation of state and the chemical 
potential remain unchanged with respect to the standard case. On one hand, the Helmholtz free 
energy does not depend on the length of the system, and, on the other, its dependence on M is 
not modified with respect to the standard case. 

The expressions for the entropy, internal energy and heat capacity are the following: 



Mk 



In I Yl [e~^ an + e ~^ K ^ 



+ 



fifrbj e "a 2 " 

~z W~ 



U 



_ gnu 

Afhui e 
~Z A 2 " 



(3.9) 



C v _ e 
Mk ~ 



Z 



+ B n e 



A 4 



E 



+1 [B n 



2A 2 



e ^2 



Z 



2A 2 

f3Hu 



2A 2 

f3hu 



n ) 



(3.10) 



where 



A,, 



-a n + 1 



B„ = ( A , b n . : + 1 j • 



These sums can not be reduced to a simple form but they can be treated numerically. These 
thermodynamical functions for arbitrary A differ significatively with respect to the standard 
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yVk 



kT 

ftco 



kT 
tuj 



(a) 



(b) 



Figure 2. Here we plot the internal energy U and the heat capacity as functions of kT/hui. In both 
graphics the solid (blue) line corresponds to the standard case with A = 0, the dashed (green) lines 
corresponds to (3.9) and (3.10), respectively, for A = 0.2. The dotdashed (red) lines correspond to the 
approximated quantities (3.14) and (3.15) also for A = 0.2. For the internal energy U we can see that the 
ground energy is shifted in the inbox of (a) to a lower value, for high temperatures the exact modified 
behavior decreases, while the approximated one increases with respect to the standard case. For the 
heat capacity Cy we notice that for low temperature the three cases are very similar (b), while for high 
temperatures the behavior is completely different. 

case (see Fig. 2). However, for small A the polymer and standard cases behave qualitatively in 
a similar manner. 

Next we study the afore mentioned limiting cases. The case A <C 1, corresponds to small 
deviations from the standard case A = 0. Although the case A > 1 has no direct interpretation 
in terms of length scales, we can say that it corresponds to the case when the energy of the 
system is much greater than the energy associated with the polymer scale. 



3.1.2 Approximate spectrum A <C 1: ensemble of harmonic oscillators 

Using the approximation given in (2.10) for small A, consistent with n max ~ A -2 , the partition 
function becomes 



A' 2 

£' 

n=0 



-/8fiw(|+n) 



cxp 



16 



n(n + 1) 



summations can be performed in a closed form to yield 



e 2 



1 



-/3fiw 



A 2 

1 + — (3hu , 



(\ 



(3.11) 



We recognize the first term in (3.11) as the usual partition function of the standard oscillator; 

_ i 

the second term is the leading polymer correction and we are neglecting terms of order 0(e a 2 ) 
which clearly tend to zero in the limit A — > 0. 

The regime where the first term is much larger than corrections, i.e. the classical regime, 
corresponds to /3Hlj S> 1. This can be seen in Fig. 2. The polymer oscillator approximates the 
standard oscillator for low temperatures. 

Using the partition function (3.11) we calculate the corresponding thermodynamic quantities, 
from (3.2). Then we can write a modified Helmholtz free energy 

2^ 



j3hu 



+ ln(l - e 
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^ + ln(l-e-^)-A^f 1 + e "2V" 
2 v ; 32 yi-e-PtoJ 



(3.12) 



Substituting (3.12) into (3.4), (3.5) we obtain: /x = F/M and p = 0: i.e. the equation of state 
does not change. As for the entropy we have 



S = Mk 



f3hw 



ln(l 



1)3 



For the internal energy we found 



U = Mhuj 



1 1 

2 ~*~ e /3fiaj . 



A 2 (l + e^)(e : 
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and for the heat capacity 



C v = Mk(j3huj)' 2 



f3tlh) 



r e 0hu — 1)2 
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A 2 (2 + /3M1 + 4e/3Sw ) + e 2l3hw {/3Hio - 2) 
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(3.13) 



(3.14) 



(3.15) 



We can see from (3.15) that heat capacity is increased due to the polymer correction, while it 
decreases the energy (3.14). Since fihuj 3> 1, such modifications are very small. 

Physically, an ensemble of harmonic oscillators can be used to model vibrations in solids 
(phonons). In a solid the vibrational modes are modeled as a collection of harmonic oscillators 
in the so-called harmonic approximation which consists of approximating the classical interac- 
tion Hamiltonian of the atoms in a solid by a second order Taylor series [59]. The simplest 
model is called the Einstein model which assumes that all vibrational modes have the same 
frequency cj, and that the oscillators are independent with no interaction. The thermodynamic 
magnitudes (3.12), (3.13), (3.14) and (3.15) contain modifications to the thermodynamics of an 
Einstein solid by defining the vibrational temperature 8y = hui/k. Note that in this case the 
polymer scale is not involved in the definition of vibrational temperature. 

3.1.3 Approximate spectrum A ^> 1: ensemble of rotors 

Finally we consider the A > 1 case for which the partition function becomes 



e a2 -f e a2 



E 2 

n=l 



e * 



(3.16) 



where $3 is the Jacobi's elliptic theta function [49]. As in previous cases, chemical potential and 
pressure are unaffected, while entropy internal energy and heat capacity are as follows: 
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Figure 3. Cy for the polymer oscillator. Blue line corresponds to the standard case with A = 0, dashed 
line is the approximation using the partition function (3.16) for A > 1, and red line corresponds to the 
case using the full spectrum in terms of Mathieu characteristic functions in the same regime. We observe 
that, except for the value of the maximum of Cy, the approximation is very good for most temperatures. 



Again we omitted the arguments and e sr of 773, while the primes in v 3 and $3, correspond 

to the first and second derivatives of $3 with respect to its second argument e s . Notice from 
Fig. 3, that the heat capacity of this system tends asymptotically to its classical value at high 
temperatures, while for low temperatures tends to zero, before passing through a maximum. This 
same behavior is found in the rotational contribution of particles with internal structure or an 
ensemble of rotors [59]. Hence for large A the quantum polymer oscillator approximates the rotor. 

Actually, one can define the corresponding rotational temperature for this system. In the 
standard case the rotational temperature is r = ^fk, with I the moment of inertia. Using 
in this case the effective moment of inertia I e ff = — p-, the rotational temperature turns out 

to be B r = ^jf- = m g^ , which does depend on A. Note that the vibrational and rotational 
temperatures for the system differ by a factor Q r /Qy = A 2 /8, which indicates that for small A 
the rotational states are negligible, while the oscillatory states are for large A. 

The full behavior of the polymer system can be seen in Fig. 2. The ensemble of polymer 
oscillators behave just as an ensemble of quantum pendulums. For low temperature it approaches 
the usual ensemble of oscillators, whereas at high temperatures it approaches an ensemble of 
rotors. This behavior is evident when studying the limiting cases separately as we have done here. 

As usual in ensembles of systems from a spectrum containing n 2 , there is a maximum in Cy 
separating different behaviors. From energy considerations, this peak appears due to a change 
of concavity in U, which can be seen in Fig. 2(a). In the polymer oscillator as well as for 
the pendulum, the maximum in Cy separates the rotational from the oscillatory behavior as 
a kind of smooth phase transition [62]. Given a fixed A and for low temperatures, the oscillatory 
states are turned on. As temperature increases more energetic states appear and for the value 
of temperature at which Cy is maximum rotational states arise. For high temperatures, only 
the rotational states remain excited. Moreover, it should be noticed that the maximum of Cy 
depends on A. For A small enough, the temperature at which the maximum occur, would be 
very high so that the system approximates very well the oscillatory behavior, while to reach the 
rotational behavior one would need a huge amount of energy. 
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3.2 Partition function for the polymer ideal gas 
3.2.1 Exact spectrum 

In the previous section we found that the energy spectrum of a polymer quantum particle in a box 
is proportional to cosn7rA, namely equation (2.21). We also showed that if we expand cosnvrA, 
when A 1, i.e., if the dimensions of the box are large compared to the minimum scale [aq, we 
obtain the usual quantum spectrum with modifications of order 0(A 2 ) equation (2.21). 

We can calculate the partition function with either the exact (2.21) or approximated spec- 
trum (2.22). Let us use the full spectrum (2.21). One way of introducing the approximation 
n < A , for which the analysis is valid, is to consider it as a cut-off in the energy, similarly as 
what has been done in [21]. First we rewrite the partition function as follows 

A 2 00 / A2 \ 

Z(f3) = e 2 ^o ex P 2 cosn7rA ' ( 3 - 17 ) 

n=0 \ 27r V J 



where A = y 2n ^[ l is the thermal wave length [59], and the exponential in the argument of 
the sum can be written conveniently as an infinite sum of modified Bessel functions of first 
order [49] 



-^-cwrwrA / A 2 \ A / A 2 \ 

H^J +2 gH^J C08(fcrarA) - (3 - 18) 

Then we can replace (3.18) in (3.17) and perform the sum over n. However, we immediately 
see that the first term diverges. It is then necessary to consider the sum only up to 1/A, as 
dictated by the approximation. Moreover, when A — > 0, the limit of the sum tends to oo, as in 
the standard case. Thus, the sum in (3.17) becomes 



M 2 A U/AM J A 2 \ ^ / a 2 \ fkTT\\ . , 

> exp k cos nirA = —in n + cosh ^ + > hi o cot sin kir. 

\ 27r Mo / A °V2vr^; V2vr/i 2 ; ^ W^J \ 2 J 

In the last term cot kn\/2 ~ 2/(kirX) for small A, then it is zero for any integer k. The partition 
function for this case has the form: 

_ A 2 

_ A 2 

Let us notice that the last term e tends to zero as //q —> 0, also there is a constant term 1/2, 

which only redefines the scale of Z 4 . Consider then the partition function as 

_ A 2 

/ A 2 \ e 2nfl 'o 



4 Moreover, one can think that this term comes from the way we approximate the sum. We can realize this 
by using the Euler-Maclaurin formula to calculate the sum in the partition function, this calculation is usual in 
standard statistical mechanics 

N , N 



n=0 J ° 



/(»)** +^(/(o) + /(#)) + • 



using N = 1/A and f(n) = exp ( cosn7rA^ we find, as we shall see, that the partition function coincides with 
the asymptotic expansion of (3.19) and also contains the factor 1/2, then it is not due to polymer corrections and 
we can ignore it from now on. 
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Figure 4. The Cy for the polymer ideal gas has a maximum around T ~ 28 po i y , for smaller values of 
temperature one recovers the usual behavior, and when T — > 00, Cy tends to zero. 

and let us calculate the corresponding thermodynamical quantities by recalling that the Helm- 
holtz free energy in the case of indistinguishable particles can be determined by the relation 



In all previous expressions we omit the arguments of the Bessel functions. We notice that also 
in this case, with the assumptions that we made, the equation of state remains unchanged. 

From Fig. 4, we notice that for low temperature (with generic (io) Cy features the usual ideal 
gas behavior namely, it takes a constant value. On the other hand as temperature increases Cy 
reaches a maximum (which depends on fio) and then goes to zero as T — > 00. This is consistent 
with the asymptotic behavior of the energy in the same regime. The maximum of Cy occurs at 
about T ~ 0p O i y /2, where © po i y = k ^ 2 = jf- Remarkably, since the energy of the polymer 
particle in the box is bounded between and E*, we are just regaining the so called Schottky 
effect, [59], that appears for a two level quantum system in which the peak of Cy appears for 
a temperature given by the energy difference between levels divided by k. 

Moreover, it is well known that systems which are bounded from above, as the present case, 
allow the existence of negative temperatures [63]. From the thermodynamical definition of tem- 
perature, negative values of T correspond to negative values in the slope of the graph of energy 
versus entropy. To see how this may happen one notices that in the partition function (3.1), and 
for positive temperatures, higher energy states contribute less than the low energy ones, this 
situation gets reversed for negative temperatures and in this situation it is mandatory that the 
energy be bounded above for the partition function to make sense. That's why only systems 





the other thermodynamical quantities are as follows: 
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such as those having two levels as magnetic systems and nuclear spin systems, feature negative 
temperatures [63]. Negative temperatures correspond to energies that are in principle experi- 
mentally accessible, due to the fact that Cy remains positive for negative temperatures, e.g. in 
nuclear spin systems [63]. 

Now let us analyze our polymeric gas in the region of negative temperatures. In the limit 
T — > + , U reaches its minimum value which is zero; however, in the limit T — > 0~, U takes 
a value that is twice its asymptotic value. In the positive temperature region the energy increases 
from zero to its maximum value as T — )• oo. In the negative temperatures regime U 

decreases from twice to once its asymptotic value as T — > — oo. We can see this behavior 
from the Fig. 5. In the polymer case, when hq is very small, the maximum energy for positive 
temperature tends to infinity and we can not reach the negative temperature regime. To access 
negative temperature in the polymer case will be very difficult in the case A <C 1. 



3.2.2 Approximate spectrum 

Let us consider the limit A — > of the partition function (3.19) in order to regain the results for 
the standard ideal gas. We can make use of the asymptotic expansion of Jo [49], which yields: 

+ -(l + e ^5 J. (3.20) 

Of the two terms in round brackets in (3.20) only the 1/2 remains in the (j,q — > limit. However, 
it only yields a constant shift in Z. 

Interestingly, a different approximation from the above is to use the approximated spec- 
trum (2.22) instead of (2.21) in the partition function. Since A C 1 we consider only the first 
terms in the series of the second exponential, so that the partition function could be expressed as 

Z^Ee-S^m + A^^,V^ +0(A 4 ). (3.21) 

n=0 n=0 

Now, by the use of the Poisson resummation formula, namely 
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we can write the partition function as 
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that is precisely the first term in (3.20). The first term is the standard partition function for 
the one dimensional ideal gas. This leads us to the constraint for the approximation to hold, 

that no < Ayf . 

Furthermore we can compute it with more orders in Ho 
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(3.22) 



Let us use (3.22) to order Ho to obtain the approximate thermodynamic quantities: 

1 



l-lnAA + ln,--ln, + ln V2 ^ + ln(l + f^ 



We notice that F diverges for f3 = or T — > oo, but with the opposite sign than in the usual 
case. As we approach to (3 = we found a maximum after which the free energy diverges to — oo 
at /3 = 0. This means that the effect of the minimum length scale on the free energy, is to 
provide a turning point for F at very high temperatures. Furthermore, when we increase the 
value of no, the free energy becomes negative at high temperatures. 

We can calculate gas pressure through the relationship p = —dF/dL, which is valid for one 
dimensional system. Realizing that the term that comes from the polymer correction does not 
depend on L, but on ho, we can obtain the equation of state for the ideal gas as p = N//3L, 
which is the same as in the continuum case. The polymer corrections can not be seen as an 
effective interaction among the particles as it was the case for GUP or MDR [35]. 

Recalling (3.3) and (3.5) the chemical potential and the entropy of the gas can be calculated 
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(3.23) 



The last expression (3.23) would be the equivalent of a one dimensional Sackur-Tetrode's formula 
for the entropy with corrections due to the underlying discreteness. The energy and heat capacity 
are obtained from (3.6) and (3.7) respectively as 



U = 
C v 



Hn TU 



AT 
2/3 
kM 

~T{ L+ 2 h 2 p 



4 h 2 /3 
1 + ^ + 



(3.24) 



All the above expressions are also obtained when considering the asymptotic expansion of Bessel 
functions of the quantities obtained in the previous section. 

We notice from (3.24) that heat capacity has a different behavior with respect to the standard 
case for which Cy is constant. For (3 = 0, or high temperature, indeed it diverges in this 
approximate case. However, we know that this is only an approximation which corresponds, 
in the full case, to increasing Cy to its maximum. Because heat capacity is related to energy 
fluctuations, we can say that for high temperatures there are strong fluctuations in the energy 
due to space discreteness. 
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Notice that if we use the result (2.23) the calculations of the thermodynamical quantities 
remain the same and only differ by numerical factors, i.e. in the partition function (3.21) a fac- 
tor 1/24 is replaced by —1/3, etc. The main difference is that with (2.23) the thermodynamical 
quantities will decrease instead of increase as our results show. 

4 Discussion 

Polymer quantum mechanics considers mechanical models quantized in a similar way as loop 
quantum gravity but in which loops/graphs resembling polymers are replaced by discrete sets 
of points. It has allowed to study some features of loop quantum gravity in a simpler context, 
namely through the use of mechanical systems [19, 20]. Indeed this opened up the possibilities 
to investigate different physical problems [19, 26, 39, 27, 28]. On the other hand important 
gravitational systems like the cosmos itself and black holes necessarily involve thermodynamics 
in their description. However, little attention has been given to the thermostatistics of poly- 
mer quantum systems. In this work we embarked on this task using the canonical ensemble 
theory applied to a polymer solid and a polymer gas, both in one dimension. The resulting 
thermodynamic quantities have modifications due to the minimum length scale that introduces 
the polymer quantization. Thermodynamics with modifications due to quantum gravity has 
been studied in the context of extensions of the standard model with Lorentz violations [36] 
and also models that mimic the graphic states of loop representation, but which, however, do 
not represent physical situations [64]. It is important to stress that here the canonical partition 
function is used as the simplest case, the relation between the discrete space and the statistics 
is an open issue [60, 61]. 

First we consider an ensemble of polymer oscillators, and we noticed that it can be interpreted 
as an ensemble of quantum pendulums with two regimes: For small A one has the vibrational 
or simple oscillator regime, and for large A the rotational regime. Both behaviors are separated 
by a maximum in the Cy. In the generic A case, we observe that to access the rotational regime 
requires a high temperature that also depends on A. In that case the partition function and 
therefore the thermodynamic quantities, are written in terms of certain sums of the exponentials 
containing the Mathieu characteristics functions (3.8). 

In the vibrational regime A <C 1 the partition function can be expressed as a power series 
on the minimum polymer length scale. Interestingly the equation of state is not modified. As 
is well known, with this model one can model an Einstein's solid introducing the vibrational 
temperature Qy that takes into account the characteristic energy of the system in this regime. 
We recover the known thermodynamical quantities when A — > 0. 

The rotational regime A > 1 that corresponds to an ensemble of rotors, can be character- 
ized by the rotational temperature G> that in this case depends on A, as r = f^A 2 . In this 
case Z and the thermodynamical variables can be written in terms of the Jacobi's theta func- 
tion $3, (3.16). This regime can be of interest for the polymer analogous to the trans-planckian 
problem, since this case corresponds to energies beyond the characteristic polymer energy given 
by the coefficient of equation (2.5). 

Amusingly it is possible to consider the effect of the polymer quantization in the case of 
electromagnetic radiation as follows. Let us recall that an ensemble of oscillators can be used to 
model excitations in solids (phonons), but also quantum excitations of the electromagnetic field 
(photons). Historically, equilibrium radiation was first studied by Planck in 1900, who conside- 
red this system as an ensemble of harmonic oscillators with the same frequency. Nowadays we 
consider that photons are ultra relativistic bosons with some particular energy hui [59]. Let us 
use Planck's simple model. By interpreting the internal energy as N times the average energy 
of each oscillator, i.e. U = J\f(E n ), from which it follows that (E n ) = %^ + (n), where (n) is the 
average occupation number that can be obtained from the previous equation and (3.14). This 



Statistical Thermodynamics of Polymer Quantum Systems 



19 



yields 
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where the second term is the polymer correction. With (4.1) we can calculate the spectral 
density that is defined as 

d fU\ . . _ g(u) fl , 

where the density of states <?(w) = ^j, V being the volume. The expression (4.2) is thus the 
black body distribution now containing polymer corrections. Note that the spectral density for 
a fixed temperature is modified for high frequencies. We know that the cosmic background 
radiation CMB with T = 3K is the black body that has been measured more accurately [65]. 
This further constrains the possible value of A and therefore the minimum scale of /xo- I n this 
case such changes do not alter the functional dependence on temperature as in other approaches 
[34, 35, 37]. Using (4.1) and (4.2) we can obtain directly the energy density by integration. It 
follows that the Stefan-Bolztmann law takes the form 

Thus, given (4.3) the Stefan-Boltzmann constant is modified by a term of order 0(A 2 ) as 

^ = ^3 (1 + 0.41648^). 

where A is small. Notice that, as expected, we recover the usual formulae of thermodynamics 
of black body radiation for A — > 0. Similar analysis using other proposal have been given 
in [34, 35, 37]. 

As for the ideal polymer gas, the partition function was determined in two regimes depending 
on whether the spectrum is considered in its exact or approximate form. In the exact spectrum 
case we notice there is a cut-off in the energy levels proportional to 1/A. This ensures in 
particular the convergence of the partition function (3.17). We note that in this system is quite 
evident that negative temperatures are allowed. In the approximate spectrum case we express the 
partition function as a series in powers of /io- Of course standard thermodynamics is contained 
in our results in the limit A — > 0. These results resemble those obtained from GUP [32, 58] in the 
sense they correspond to quadratic corrections although they feature opposite tendencies [44]. 
Clearly within our analysis the origin of the modifications can be traced all the way back to the 
polymer model, in which there is also a modified uncertainty relation [39]. 

In regard to the thermodynamics of both of our systems it is worth stressing a common 
behavior of their heat capacity. For the gas, in contrast to the standard case, the Cy has 
a maximum of about half of O po iy, the temperature at which the polymer effects are evident. 
The polymer oscillator, on the other hand, also shows a similar behavior in its Cy, having 
a maximum that separates two different behaviors characterized by 0y and r . 

Now we mention some possible extensions of this work. The canonical distribution adopted 
here is only intended to be an approximation. The statistics of polymer systems, which are 
naturally discrete is still an open issue [60]. In [61] the problem of calculating the number 
of accessible microstates was considered in a semiclassical perspective; it was found that the 
standard methods yield only approximate results. To count states in polymer systems and to 
include fermions and bosons, again a better understanding of the statistics of polymer quantum 
systems is needed. 
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In this work we learnt that polymer quantization induces some effects on the statistical de- 
scription of the systems, which in principle we can explore. This effects are such as the maximum 
in Cy between two different behaviors, or the appearance of negative temperatures in the case 
of ideal gas. Interestingly, a similar behavior has been observed in gravitational systems such 
as black holes, where the heat capacity shows a phase transition that depends on the minimum 
length scale [66] . The study of the thermodynamic quantities of gravitational systems has fo- 
cused mainly on obtaining the entropy of black holes from gravitational quantum states [67], or 
to find effective modifications to those thermodynamical quantities, [41, 68]. However, recently 
there has been some quantum gravity models which are based on thermostatistics of condensed 
matter systems, can yield interesting results [69, 70, 71]. 

Certainly an interesting question emerges when exploring the thermodynamics of non equi- 
librium processes such as those that occur in the early universe or in the black hole evaporation. 
Clearly it is necessary to extend the methods presented here to address these problems. 
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